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ABSTRACT 

Spectral  and  cross-spectral  analysis  are  used  to  investigate 
the  existence  and  structure  of  wave-type  disturbances  in  the 
troposphere  over  the  equatorial  South  Atlantic  Ocean,  The 
data,  which  came  from  the  NMC's  tropical  analysis,  consists  of 
wind  components  for  5  levels  of  the  troposphere  (700  to  200  mb) 
and  comprise  a  period  of  411  days.   Also  data  at  the  stations 
Parnaiba  Airport  r(3°S,  42°W)  and  Ascension  Island  (8°S,  14 °W) 
are  examined.  The  spectral  results  indicate,  with  95$  confidence, 
that  there  is  periodicity  in  the  time-series  analysed  in  the 
frequency  band  0ol6-0,24  cpd  (approximately  4-  to  6-  day  periods). 

Cross-spectral  analysis  results  indicate  that  the  disturbance 
has  a  wavelength  of  about  6000  km,  propagates  westward  at  a 
rate  of  about  11°  longitude  per  day,  has  a  slight  NW-SE  tilt 
and  tilts  backward  in  the  vertical. 

The  energetics  of  the  disturbance  could  not  be  determined. 
Examination  of  the  shear  of  the  wind  field  and  the  tilt  of  the 
disturbance  indicates  that  barotropic  instability  is  unlikely 
to  be  an  energy  source. 
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I.   INTRODUCTION 

Wave-type  disturbances  of  synoptic  scale  in  the  tropical 
regions  have  been  extensively  studied  in  recent  years.   Among 
others,  the  works  of  Yanai  et  al.  (1968),  Wallace  and  Chang 
(1969),  Chang  (1970a),  Chang  et  al.  (1970b),  Nitta  (1970),  Reed 
and  Recker  (1971)  and  Wallace  (1971)  for  the  Central  and  Western 
Pacific  and  Burpee  (1972,  1974)  for  North  Africa  have  documented 
the  existence  of  such  disturbances  with  periods  ranging  from 
four  to  six  days  and  wavelengths  from  2000  to  8000  km.  These 
studies  used  either  the  spectral  analysis,  composite  techniques 
or  a  combination  of  both  techniques. 

Due  to  the  absence  of  an  adequate  observational  network, 
the  weather  over  the  equatorial  South  Atlantic  Ocean  is  generally 
represented  by  its  mean,  undisturbed  state.  There  are  few 
references  in  the  literature  about  the  perturbations  caused  by 
organized  synoptic  systems.   Yet,  there  is  evidence  of  dis- 
turbances crossing  the  west  coast  of  Africa  into  the  ocean 
between  the  equator  and  5°  S  (Mbele  Mbong  1974)  and- the  east 
coast  of  South  America  into  the  continent  at  approximately  10  S 
(Ramos  1974) . 

The  availability  of  wind  data  for  five  levels  of  the  tro- 
posphere from  analysis  produced  by  the  Tropical  Analysis 
Group  at  the  National  Meteorological  Center  suggested  an  in- 
vestigation about  the  presence  of  wave-like  disturbances  in 
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the  equatorial  South  Atlantic  area.  This  tropical  analysis 
is  performed  at  12  hour  intervals  based  on  all  available  con- 
ventional radiosonde  observations,  satellite  derived  winds  and 
aircraft  reports.   Due  to  the  scarcity  of  reporting  stations, 
the  data  in  the  equatorial  South  Atlantic  are  fairly  limited. 
However,  the  NMC  analysis  may  still  contain  some  information 
of  the  synoptic-scale  disturbances  passing  through  this  area. 
Since  the  amplitude  of  the  disturbances  that  can  be  picked 
up  by  the  objective  analysis  scheme  are  expected  to  be  small 
and  easily  obscured  by  the  persistence  which  forms  a  substantial 
portion  of  the  total  variance  of  the  wind  field,  these  dis- 
turbances may  be  very  difficult  to  discern  by  analysis  of  indi- 
vidual weather  maps  on  a  routine  basis.   On  the  other  hand, 
the  chance  of  detecting  them  should  be  considerably  improved 
when  long-period  time  series  are  analyzed. 

The  purpose  of  this  study  is  to  perform  such  an  analysis 
of  the  NMC  data  using  spectral  and  cross-spectral  techniques 
in  the  hope  that  some  information  on  wave-type  disturbances 
in  the  equatorial  South  Atlantic  contained  in  the  data  can  be 
obtained. 


II.   DATA 

A.   BASIC  DATA  SET 

In  addition  to  the  NMC  analysis  mentioned  earlier,  two 
other  sets  of  data  were  made  available  during  the  course  of 
the  work.   A  description  of  each  set  follows. 

The  NMC  analysis  contains  the  zonal  (u)  and  meridional 
(v)  wind  components  for  five  levels  of  the  troposphere,  namely, 
700,  500,  300,  250  and  200  mb.   These  winds  are  available  at 
rectangular  grid  points  evenly  spaced  by  five  degrees  of 
longitude  on  a  Mercator  projection.   The  grid  distance  is 
approximately  556  km  at  the  equator  (Bedient  and  Irwin  1970). 
The  area  selected  for  this  study  is  bounded  by  45°W,  15°E, 
10°S  and  the  equator,  including  39  grid  points  with  13  at  each 
latitude  (Figure  1) *   In  total,  390  time  series  were  generated, 
each  series  containing  twice-daily  data  for  a  period  of  441 
days  (8£2  data  points)  from  1200  GMT,  23  January  1971  to  0000 
GMT,  S  April  1972. 

The  second  set  of  data,  provided  by  the  National  Climatic 
Center,  consisted  of  surface  and  radiosonde  observations  for 
a  few  stations  in  the  African  Coast  and  Ascension  Island  (3  S, 
14 °W) •   However,  only  data  at  Ascension  Island  were  analyzed 
because  the  excessive  amount  of  missing  data  at  the  other 
stations  rendered  their  time  series  unsuitable  for  spectral 
analysis.   For  Ascension  Island,  23  time  series  of  daily 
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observations  of  horizontal  wind  components,  temperature  and 
mixing-ratio  at  eight  pressure  levels  (1000,  $50,  700,  500, 
400,  300,  250  and  200  mb)  were  generated,  spanning  a  period 
of  36l  days  from  4  January  1971  to  31  December  1971  (no  mixing 
ratio  for  the  four  upper  levels).  Since  on  the  average  about 
40^  of  the  data  are  missing  in  each  time  series,  even  these 
data  were  considered  marginally  usable  for  a  spectral  analysis,, 

The  third  set  came  from  the  Brazilian  space  research 
institute  (instituto  de  Pesquisas  Espaciais)  through  the 
hydrographic  agency  of  the  Brazilian  Navy  (Diretoria  de 
Hidrografia  e  Navegacao).   It  consists  of  daily  pilot-balloon 
observations  at  height  levels  1$0  m  apart  at  several  stations 
and  radiosonde  observations  of  temperature  and  relative  humidity 
at  several  others,  all  in  the  Brazilian  coast.   However,  only 
the  lower  levels  at  the  station  Parnaiba  Airport  (3  S,  42  W) 
were  considered  to  have  sufficient  data  in  the  time  domain  to 
make  spectral  analysis  worthwhile.  Time  series  of  u  and  v 
components  were  thus  generated  at  this  station  for  each  of  the 
12  lowest  height  levels,  resulting  in  a  total  of  24  series 
encompassing  the  period  from  18  June  1971  to  31  December  1971, 
These  series  represent  approximately  the  lowest  2  km  of  the 
atmosphere, 

B.   PREPARATION  OF  THE  DATA 

After  the  time  series  were  generated,  all  missing  data 
were  linearly  interpolatedo   The  fraction  of  the  interpolation 
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for  the  NMC  data  series  is  always  less  than  8.0$,  For  Ascension 
Island  this  amount  is  increased  to  about  40$  of  the  series, 
with  missing  data  more  or  less  evenly  scattered  in  the  entire 
period.  The  amount  of  interpolation  at  Parnaiba  Airport  is 
close  to  24$  with  no  large  gaps.   These  percentages  suggest  that 
the  spectral  results  obtained  from  the  station  data,  especially 
those  of  Ascension  Island,  are  less  reliable  than  the  grid- 
point  results. 

Each  time  series  was  then  subjected  to  a  high-pass  filter, 
with  the  primary  purpose  of  avoiding  leakage  of  variance  from 
the  very  low  frequencies,  which  is  quite  predominant  due  to 
the  persistence  in  the  data,  especially  the  zonal  wind  com- 
ponent.  The  weights  of  the  filter  were  obtained  from  a  normal 
curve  with  standard  deviation  50/3  half -days  for  the  twice- 
daily  series  and  25/3  days  for  the  daily  series  (Holloway 
195$) •   The  Nyquist  frequency  is  1  cycle  per  day  (cpd)  for  the 
twice-daily  series  and  2  cpd  for  the  others.   No  measure  was 
taken  to  avoid  "aliasing"  from  frequencies  higher  than  the 
Nyquist  frequency  since  not  much  variance  is  thought  to  exist 
in  them. 

The  theoretical  response  function  of  the  filter  is  greater 
than  95$  for  frequencies  higher  than  0.05  cpd  (period  of  20 
days),  and,  for  all  practical  purposes,  100$  for  frequencies 
higher  than  0.0S  cpd  (period  of  12c5  days) o   An  estimate  of 
the  actual  response  function  was  obtained  by  determining  the 
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variance  spectrum  of  the  series  before  filtering  (G,  (f))  and 
after  filtering  (G  (f))  and  taking  the  positive  square  root 
of  their  ratio  (H  (f)=G  /G  ).   This  estimate  is  displayed  in 
Figure  2. 

After  filtering,  the  time  series  at  the  grid  points  were 
reduced  from  832  to  7$4  points*  The  Ascension  Island  series 
were  reduced  to  313  points  and  the  Parnaiba  Airport  series  to 
149  points. 

The  variance  spectrum  for  each  series  was  then  computed 
with  a  maximum  lag  of  25  days,  which  implies  50  lags  of  one- 
half  day  for  the  grid-point  series  and  25  lags  for  the  station 
series.   The  cross-spectra  were  computed  whenever  possible, 
between  series  at  different  levels  (inter-z),  different  longi- 
tudes (inter-x),  different  latitudes  (inter-y)  and  between 
the  two  wind  components. 

All  the  spectral  computations  were  performed  at  the  W.  R. 
Church  Computer  Center  of  the  Naval  Postgraduate  School  using 
the  program  BMD02T  of  the  UCLA  Biomedical  Statistical  Program 
Package. 
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III.   TEST  OF  SPECTRAL  SIGNIFICANCE 

The  significance  of  the  spectral  results  was  tested  using 
a  procedure  recommended  by  Mitchell,  et  alo  (1966).  This 
procedure  is  briefly  described  below  along  with  pertinent 
considerations • 

The  BMD02T  program  used  to  calculate  the  variance  spectrum 
begins  by  computing  the  mean  of  the  data  and  subtracting  it 
from  each  value  of  the  series.  The  serial  autocorrelation 
coefficients  (Rn)  for  lags  from  zero  to  M  are  then  calculated, 
where  M  (maximum  number  of  lags  chosen)  Z.  N  (number  of  points 
in  the  series) •  The  cosine  transforms  of  these  M+l  lag  auto- 
correlation values  are  then  computed,  yielding  M+l  "raw"  esti- 
mates of  the  variance  spectrum,  each  value  being  a  rough  measure 
of  the  total  variance  in  the  original  series  that  is  contributed 
by  wavelengths  near  each  harmonic  of  the  fundamental  wavelength 
of  the  analysis  (2M).  The  raw  estimates  are  then  smoothed  by 
a  "hamming"  window  function  (a  3  term  weighted  moving  average 
with  weights  equal  to  0,23,  0.5A-  and  oe23,  respectively),  with 
the  purpose  of  deriving  a  consistent  estimate  of  the  final 
spectrum  in  terms  of  M+l  discrete  estimates. 

To  evaluate  the  cross-spectrum,  the  program  first  calculates 
the  serial  cross-correlation  coefficients  in  the  same  form  as 
for  the  autocorrelation  coefficients.   After  that,  the  co- 
spectrum  is  obtained  by  determining  the  cosine  transforms  of 
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these  cross-correlation  coefficients  and  the  quadrature  spectrum 
is  obtained  from  the  sine  transforms  of  the  cross-correlation 
coefficients.   Both  the  co-spectrum  and  the  quadrature  spectrum 
are  smoothed  by  the  "hamming"  method.  From  the  co-spectrum  and 
the  quadrature  spectrum,  the  amplitude  and  phase  of  the  cross- 
spectrum  are  computed  and  the  coherence  square  between  the  two 
series  is  determined. 

In  this  study,  each  computed  variance  spectrum  was  tested 
by  fitting  a  "null"  hypothesis  continuum  spectrum.  The  choice 
of  a  "null"  continuum  was  based  on  the  observed  fact  that  the 
lag-one  serial  autocorrelation  coefficient  R.,  of  the  series 
analyzed  differs  markedly  from  zero,  being  typically  of  order 
of  0.70  for  the  u  series  and  0.55  for  the  v  series.   This  indi- 
cates the  possible  large  influence  of  simple  persistence  to 
the  series.   If  a  series  is  made  up  entirely  of  random  varia- 
tions and  persistence,  it  may  be  written  as 

Xi+1  =  HjXj  +£  .+1 

where  L.    -,  is  the  random  part  of  X-  -,  independent  of  any  of 
the  preceding  terms  of  the  series  and  R-iX-  is  the  part  of  X-  -, 
"explained"  by  X- .   Applying  the  same  reasoning  for  X.  ^i    e"tc> 
and  combining  the  equations  thus  obtained,  it  can  be  seen  that 
the  sample  serial  autocorrelation  coefficients  must  approxi- 
mately  bear  an  exponential  relationship,  R  =  R-,  •   It  follows 
that  the  autocorrelation  function  must  decrease  exponentially 
and  approach  zero  asimptotically  (see,  for  example,  Brooks  and 
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Carruthers  1953,  or  Mitchell  1963).   In  such  case,  the  appropriate 
"null"  continuum  should  also  be  a  decreasing  function  whose 
shape  depends  on  the  value  of  the  lag-one  autocorrelation 
coefficient.   Assuming  that  R*  is  an  unbiased  estimate  of  the 
unknown  lag-one  autocorrelation  coefficient  for  the  population, 
this  "null"  continuum  may  be  described  by 

SK  =  S  (  (1-R12)/(1+R12  -  2R1  COS  (K/M)  )  )    , 

where  3  is  the  average  of  all  M+l  "raw"  spectral  estimates  in 
the  computed  spectrum,  K  is  the  harmonic  number  and  SK  is  the 
variance  per  unit  frequency  interval. 

The  first  step  of  the  testing  procedure  is  to  verify 
whether  the  autocorrelation  coefficients  for  the  various  lags 
of  a  given  time  series  drop  exponentially  or  not.   Only  in  the 
affirmative  case,  the  chosen  model  is  considered  appropriate 
and  the  variance  spectrum  may  be  compared  to  the  "null"  continuum 
in  order  to  determine  if  the  series  contain  some  other  form  of 
non-randomness  in  addition  to  persistence.   The  value  of  the 
spectral  estimate  for  each  frequency  interval  is  compared  with 
the  local  value  of  the  "null"  continuum.   If  it  differs  from 
the  continuum  by  a  statistically  significant  amount,  the  hypothesis 
of  non-randomness  other  than  persistence  in  the  frequency  band 
examined  is  accepted. 

The  statistic  associated  with  each  spectral  estimate  (i.e., 
with  the  variance  per  interval  of  frequency)  is  the  ratio  of 
the  magnitude  of  each  spectral  estimate  to  the  local  magnitude 
of  the  continuum.   This  statistic  has  been  shown 
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2 

to  be  distributed  according  to  a  chi-square  (  %    )    divided  by 

the  degrees  of  freedom  (D0F)o  The  degrees  of  freedom  are  given 
by 

DOF  =  (2N-M/2)/1VI       , 
where  N  is  the  sample  size  and  M  is  the  maximum  lag.   It  follows 
that,  for  a  chosen  level  of  significance  (0,05  in  the  present 
study) ,  confidence  limits  can  be  established  by  multiplying 
the  "null"  continuum  by  the  proper  value  of  %  /D0Fo  The 
"null"  hypothesis  may  be  then  phrased  as  follows:   The  only 
form  of  non-randomness  in  the  series  is  persistence.   If  the 
estimates  of  variance  per  interval  of  frequency  exceed  the  con- 
fidence limits,  the  hypothesis  is  rejected  and  it  is  concluded 
that  there  exists  periodicity  in  the  frequency  band  under  con- 
sideration.  Figures  3(a)  and  (b)  illustrate  the  analysis 
procedure.   They  were  generated  from  the  1000  rnb  v-series  at 
Ascension  Island,   Figure  3(a)  shows  the  autocorrelation 

function  R   (full  line).  The  dashed  line  represents  the  expo- 

n        n 
nentially  decreasing  function  R.,  •   As  R,  ^  R  for  the  first 

few  lags,  the  series  is  considered  to  fit  the  first-order  auto- 
regressive  model,  i,e,,  to  be  affected  only  by  simple  persistence 
R-,  is  then  used  to  generate  the  "null"  spectrum,  which,  multi- 
plied  by  the  95th  percentage  point  of  the  %   /DOF,  gives  the 
confidence  limit  represented  by  the  dashed  line  in  Figure  3(b), 
Since  the  estimate  of  the  variance  at  the  frequency  0,20  cpd 
exceeds  the  confidence  limit,  it  is  concluded  that  there  is 
activity  in  the  series  at  this  frequency,  with  95$  confidence. 
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It  should  also  be  noticed  that  the  described  procedure 
applies  to  an  a  priori  chosen  frequency  interval,  which  is  based 
on  considerations  independently  derived.   In  the  present  study, 
the  frequency  band  from  0.16  to  0,24  was  selected  based  on 
the  results  from  previous  studies  in  several  other  areas  that 
used  not  only  spectral  analysis  but  also  composite  techniques 
(as  mentioned  in  the  Introduction),   The  need  for  an  a  priori 
selection  of  a  frequency  band  for  investigation  is  due  to  the 
fact  that  for  a  sample  variance  spectrum  with,  say,  100  spectral 
estimates,  on  the  average  5  of  these  values  are  expected  to 
lie  beyond  the  95$  confidence  limit.   So,  if  a  particular 
spectral  estimate  is  chosen  for  study  merely  because  it  is  the 
largest  in  the  spectrum,  the  usual  test  of  significance  is 
not  applicable,  and  a  different  principle  of  decision  theory 
has  to  be  applied  (Brooks  and  Carruthers  1953,  Mitchell  1963)* 

The  grid-point  series  have  7$4  data  points  after  filtering; 
the  number  of  lags  chosen  was  50  (25  days).  These  values  led 
to  spectral  estimates  with  approximately  31  degrees  of  freedom 
and  %  /DOF  =  I.46.  The  latter  value  becomes  1,52  for  the 
Ascension  series  (25  lags  used)  and  1,77  for  the  Parhaiba 
series. 

Out  of  approximately  400  series  analyzed,  only  less  than 
1.0$  do  not  fit  the  linear  autoregressive  "null"  continuum. 
Thus,  more  than  99$  of  the  series  were  tested  according  to 
this  procedure. 
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Confidence  limits  for  the  coherence  square  estimates  were 
also  established  using  the  probability  points  of  the  distribu- 
tion of  squared  coherence  compiled  by  Amos  and  Koopmans  (1963). 
The  degrees  of  freedom  of  these  estimates  are  given  by  1.25  N/M 
which  is  appropriate  for  the  hamming  smoothing  function. 

It  should  be  noticed  that  this  test  for  significance  of 
the  coherence  square  implies  a  "null"  hypothesis  of  zero  popu- 
lation coherence.   Such  a  model  is  thought  to  be  a  fairly  good 
assumption  in  the  case  of  the  700-  and  500-  mb  series  of  the 
grid-point  data  and  the  series  for  the  station  data,  as  they 
generally  show  very  low  values  of  background  coherence.  For 
the  upper  levels,  a  background  of  high  values  of  coherence 
were  found.  The  appropriate  confidence  limits  for  the  coherence 
square  should  be  raised  to  take  into  account  the  already  high 
level  of  background  coherence.  This  was  not  done  in  the  present 
study  because  the  variance  spectra  for  those  levels  generally 
do  not  show  significant  peaks  in  the  frequency  interval  studied. 
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IV.   RESULTS 

A.   VARIANCE  SPECTRA 

1.   700  mb  and  500  mb  Grid-Point  Series 

The  variance  spectra  of  u  and  v  at  the  700  mb  level 
for  all  grid-points  are  shown  in  Figure  4.   At  this  level, 
45  of  the  73  series  show  significant  spectral  peaks  in  the 
frequency  band  of  0,16-0 ,24  cpd,  with  24  peaks  in  the  v  com- 
ponent and  21  in  the  u  component.   In  general,  they  are  more 
pronounced  away  from  the  equator.   At  the  500  mb  level  (Figure 
5),  50  series  show  significant  peaks  in  this  frequency  band 
and  34  of  them  are  in  the  v  component.   Generally,  the  spectral 
peaks  seem  to  be  somewhat  more  prominent  as  compared  to  those 
at  700  mb.   In  addition,  no  strong  latitudinal  preference  is 
detected  at  this  level. 

In  order  to  better  establish  the  significance  of  the 
spectral  peaks  for  the  0,16-0,24  cpd  frequency  band,  a  composite 
spectrum  was  constructed  for  each  wind  component  and  level  by 
averaging  the  spectra  at  the  13  grid-points  of  each  latitude. 
The  composite  spectrum  was  then  compared  to  an  "average" 
"null"  continuum  calculated  from  the  average  of  the  lag-one 
correlation  coefficient  of  the  13  series.   The  resultant  curves 
are  shown  in  Figure  7,   It  can  be  readily  seen  that  the  v  com- 
ponent has  a  well  defined,  significant  spectral  peak  in  the 
frequency  band  being  considered  at  500  mb  at  all  three  latitudes. 
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At  700  mb  the  spectral  peaks  are  less  pronounced  but  significant 
variance  in  the  v  component  still  shows  up  at  5°S  and  10°S  for 
the  0,20-0,24  cpd  frequency  band.  Thus  these  results  clearly 
indicate  the  significance  of  oscillations  of  the  v  component 
in  the  frequency  range  of  0,16-0,24  cpd,  which  corresponds  to 
a  quasi-periodicity  of  approximately  4-6  days.   On  the  other 
hand,  no  corresponding  peaks  in  the  composite  u-spectra  can 
be  found,  except  for  a  slight  indication  of  variance  at  500 
mb  level  at  the  equator.  Therefore  a  major  portion  of  the 
variance  of  the  4-6  days  oscillations,  at  least  as  reflected 
in  the  NMC  analysis,  seems  to  be  confined  to  the  v  component. 

The  frequency  band  0,16-0,24  cpd  is  the  a  priori 
chosen  band  of  interest  for  this  study.   In  the  subsequent 
sections  all  the  cross-spectral  results  discussed  are  the  average 
values  for  this  frequency  band, 

2,   Upper-Level  Grid-Point  Series 

The  number  of  significant  spectral  peaks  at  the  three 
upper  levels  (200,  250  and  300  mb)  of  the  NMC  analysis  is 
greatly  reduced.   Furthermore,  the  spectra  for  both  wind  com- 
ponents at  these  levels  resemble  each  other,  having  the 
characteristics  of  a  red-noise  spectrum  because  most  of  their 
variance  is  in  the  very  low  frequencies.  The  lag-one  auto- 
correlation  coefficients  for  these  series  range  between  0,65 
and  0,75,  somewhat  higher  than  those  for  the  series  at  500 
and  700  mb.  Thus  persistence  seems  to  dominate  the  total 
variance  at  the  upper  levels.   Since  inter-level  cross-spectra 
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results  (to  be  discussed  later)  also  indicate  a  very  high 
coherence  between  these  levels  at  all  frequencies,  it  is  felt 
that  they  are  mainly  constructed  from  extrapolations  of  a 
single  data  level.  For  this  reason,  only  the  results  of  the 
200  mb  level  are  shown  to  represent  the  upper  troposphere. 

The  variance  spectra  at  200  mb  for  each  grid-point  are 
depicted  in  Figure  6,   It  is  clear  that,  except  for  a  few 
isolated  grid-points,  mainly  at  10  S,  no  significant  spectral 
peak  can  be  found.  The  composite  spectra  at  this  level,  which 
is  also  included  in  Figure  7,  show  no  significant  activity  at 
all.   It  is  therefore  concluded  that  the  4-6  day  oscillations 
are  mainly  confined  to  the  middle  and  lower  levels  of  the  NMC 
data, 

3,   Station  Series 

At  Parnaiba  Airport,  spectra  of  u  and  v  series  for  the 
12  height  levels  from  surface  to  2070  m  were  computed  from  the 
available  pilot-balloon  data  (not  shown).  The  v  series  have 
significant  spectral  peaks  in  the  0,16-0,24  frequency  band  in 
10  of  the  12  levels,  while  the  u  series  exhibit  peaks  only  at 
the  4  lowest  levels  (surface  to  612  m) ,   Because  of  the  strong 
similarity  of  the  spectra  at  the  various  levels,  a  composite 
spectrum  for  v  series  was  generated  by  averaging  the  12  v 
spectra  in  the  vertical  and  compared  to  a  "null"  continuum, 
in  the  same  manner  earlier  described  for  the  composite  grid-point 
spectra  (Figure  8),   A  distinct  spectral  peak  centered  at  the 
0,20  cpd  frequency,  or  the  5  day  period,  is  clearly  visible. 
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For  Ascension  Island,  spectra  of  the  u  and  v  components, 
temperature  and  mixing-ratio  at  the  eight  previously  mentioned 
pressure  levels  were  computed.  The  spectra  of  the  u  series  do 
not  show  any  significant  peak  and  only  the  1000- mb  v  spectrum 
has  a  significant  peak  around  the  0.20  cpd  frequency  band 
(Figure  3)»   The  temperature  and  mixing-ratio  series  exhibit 
very  small  fluctuations  and  their  spectra  (not  shown)  resemble 
that  of  a  random  series,  with  no  feature  being  distinguishable 
above  the  background  noise. 

The  variance  spectra  for  Ascension  Island  are  expected 
to  resemble  those  for  the  neighbouring  grid  points.   Since  this 
is  not  quite  the  case,  a  careful  comparison  between  the  data 
at  Ascension  Island  and  those  at  the  two  grid  points  just  north 
and  south  of  the  island  (5°S,  15°W  and  10°S ,  15°W)  was  performed. 
A  visual  inspection  of  the  unfiltered  700  mb  u  and  v  series  at 
the  grid  points  for  both  0000  and  1200  GMT  show  them  to  be  in 
fairly  good  agreement  with  the  once-daily  series  at  Ascension 
Island,  except  that  the  latter  series  are  much  smoother  due  to 
the  excessive  interpolation  of  missing  data.  The  higher  frequencies 
at  Ascension  Island  are  therefore  strongly  dampedo   Cross-spectra 
between  the  Ascension  Island  and  grid  point  series  were  also 
computed.  The  coherence-square  values  are  generally  large  for 
the  lower  frequencies  and  decrease  towards  higher  frequencies 
but  remain  above  the  95$  confidence  limit  for  frequencies  around 
0.20  cpd.   These  facts  support  the  hypothesized  cause  for  the 
apparent  disagreement  between  the  grid  points  and  the  station 
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results,  namely,  the  excessive  amount  of  missing  data  at 
Ascension  Island  (approximately  40$). 

B.   CROSS-SPECTRA 

1,   Inter-Level  Cross-Spectra 

The  inter-level  (inter-z)  cross-spectra  for  the  grid- 
point  wind  components  were  determined  using  the  500  and  250 
mb  levels  as  the  base  series.   The  cross-spectra  with  250  mb 
as  the  base  series  reveal  very  high  values  of  coherence  square 
at  all  frequencies  between  this  level  and  the  two  other  upper 
levels  (200  and  300  mb)c  The  phase  differences  are  also  very 
small  across  the  spectrum  at  all  grid  points0  This  result, 
plus  the  previously  mentioned  resemblance  of  the  variance 
spectra  for  these  levels,  are  attributed  to  the  special  pro- 
cedure used  by  the  NMC  analysis.   In  this  region,  the  data 
are  usually  available  at  only  one  upper  level  at  each  location, 
which  may  be  from  either  aircraft  reports  or  satellite-derived 
winds.  These  data  are  then  extrapolated  to  the  other  levels, 
resulting  in  the  strong  similarity  between  the  time-series  at 
the  three  upper  levels  at  most  grid  points. 

The  cross-spectra  between  the  200,  500  and  700  mb 
series  for  both  wind  components  are  shown  in  Figure  9*   At  each 
grid  point  a  line  is  used  to  indicate  the  vertical  tilt  of  the 
disturbance  axis  between  the  700,  500  and  200  mb  levels.  This 
line  was  drawn  according  to  the  phase  differences  between  the 
500  mb  level  and  the  other  levels.   Only  phase  values 
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associated  with  a  coherence-square  value  that  meets  or  exceeds 
the  95^  confidence  limit  are  plotted.   It  is  seen  that  only 
about  1/3  of  the  grid  points  show  significant  inter-level  cross- 
spectral  values.   Surprisingly,  the  relation  between  the  500 
mb  and  200  mb  is  somewhat  better  than  that  between  500  mb  and 
700  mb.   For  those  phase  values  associated  with  significant 
coherence  squares,  the  lower  level  usually  leads  the  upper 
level,  suggesting  an  upward  propagation  direction.  This  is 
generally  true  for  both  components,  except  for  the  southeastern 
corner  of  the  region  where  the  v  component  shows  a  downward 
propagation  direction  between  500  and  700  mb.   In  almost  all 
cases  the  phase  differences  between  two  "adjacent"  levels  are 
less  than  0.20  cycle. 

The  inter-z  cross-spectra  for  Parnaiba  Airport  between 
the  12  levels  from  surface  to  2070  m  (  B00  mb)  were  computed 
with  the  SOl-m  level  series  as  the  base  series0   The  resultant 
coherence  squares  for  u  are  generally  below  the  95/£  confidence 
limit.  The  coherence  squares  for  v,  on  the  other  hand,  all 
exceed  this  limit  and  the  phase  differences  are  plotted  in 
Figure  10.   It  can  be  seen  that  the  waves  propagate' upward 
above  the  1170  m  level  (  3&0  mb) ,  consistent  with  the  results 
obtained  for  the  grid-point  series.   Below  this  level,  the 
waves  show  a  downward  propagation  direction  with  some  variations 
between  801  and  414  m.   In  general  the  phase  differences  are 
quite  small  being  less  than  0.20  cycle  between  all  levels 
except  the  surface. 


2.   Inter-Longitude  Cross-Spectra 

At  each  level  the  time-series  of  each  component  at 
15°W  were  used  as  the  base  series  to  cross  with  the  other  grid- 
points  at  the  same  latitude.   The  resultant  inter- x  cross- 
spectra  are  plotted  in  Figure  11.   In  each  panel  the  phase  dif- 
ference with  coherence  value  satisfying  the  criterion  of  95$ 
confidence  are  plotted  as  a  function  of  longitudinal  distance 
between  the  base  series  and  the  other  serieso  With  only  a 
few  exceptions,  usually  occuring  at  the  easternmost  grid- 
points,  the  coherence-square  values  between  longitudes  are 
generally  quite  high.   At  all  latitudes  and  levels  the  v  com- 
ponent shows  a  quite  consistent  east-west  phase  relationship. 
The  eastern  series  always  lead  the  western  ones,  indicating  a 
westward  propagation  direction.   The  slope  of  the  plotted 
curves  are  also  very  similar  and  show  a  longitudinal  difference 
of  about  55°  or  approximately  6000  km  for  a  complete  cycle • 
This  value  thus  represents  the  east-west  wavelength  of  the 
waves.   The  u  component  generally  indicates  a  smaller  phase 
difference  between  longitudes  and  not  so  consistent  from  level 
to  level  as  it  is  for  the  v-series,  although  the  lower  levels 
still  indicate  a  westward  propagation  direction.   Since  the 
v  component  has  most  of  the  variance  in  this  frequency  band, 
the  result  of  6000  km  is  taken  as  the  estimate  for  the  horizontal 
wavelength. 
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3.   Inter-Latitude  Cross-Spectra 

The  inter-latitude  (inter-y)  cross-spectra  for  both 
components  are  presented  graphically  in  Figure  12„   Here  each 
line  connecting  three  grid  points  at  the  same  longitude  and 
level  indicates  the  phase  differences  between  5°S  (the  base 
series),  the  equator  and  10  S.  Most  u  series  and  all  v  series 
have  significant  inter-y  coherence  squares  that  meet  the  95$ 
confidence  limit .  With  the  exception  of  a  few  cases,  the  u 
series  tilt  NW-SE,  the  average  phase  difference  between  the 
equator  and  10°S  being  close  to  0„lS  cycle.  The  v  series 
generally  have  very  small  phase  differences  between  latitudes 
with  the  average  phase  difference  indicating  a  very  slight  NW-SE 
tilt  for  the  waves.   Since  most  of  the  variance  is  in  the  v 
component,  these  phase  relationships  indicate  a  slight  NW-SE 
tilt  of  the  axis  of  the  waves, 

A-»   Inter-Parameter  Cross-Spectra 

The  cross-spectra  between  the  u  and  v  series  at  each 
grid  point  and  level  (not  shown)  have  only  about  15$  of  the 
pairs  with  significant  coherence.  For  these  pairs,  a  1/2  cycle 
out-of-phase  relationship  between  the  two  components  is 
generally  indicated.   For  non-divergent  motion  such  phase 
implies  a  NW-SE  tilt  of  the  wave  axis  (Gruber  1974). 
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V  .   SUMMARY  AND  CONCLUSION 

i  investigation  about  the  existence  of  synoptic-scale 
;ype  disturbances  over  the  tropical  South  Atlantic  Ocean, 
:he  equator  to  10°S  was  carried  out  using  the  spectral 

is  technique, 
le  main  data  base  came  from  the  NMC  tropical  analysis 

is  available  at  5°  x  5°  longitude  grid  points  and  con- 

of  winds  at  five  level's  in  the  troposphere,  spanning  a 
i  of  approximately  14  months.   In  addition,  one  year's 
from  Ascension  Island  and  a  Brazilian  coast  station, 
iba  Airport,  were  examined.  The  data  were  first  organized 
bime  series,  with  missing  data  being  linearly  interpolated, 
h-pass  filter  was  then  applied  to  the  series  to  avoid 
ge  of  the  variance  through  the  lobes  of  the  spectral 
w  in  the  computation  of  the  spectra, 

statistical  test  of  the  variance  spectra  was  performed 
suming  a  "null"  continuum  with  simple  persistence  as  the 
form  of  non-randomness.  The  a  priori  selected  frequency 
under  scrutiny  is  from  0,16  to  0,24  cpd  (4  to  6  day  periods) 

found,  with  957°  confidence  that  oscillations  are  indeed 
nt  in  the  frequency  band  in  the  meridional  component  of 
00  and  500  mb  levels  of  the  NMC  grid  point  data  and  in  the 
evel  station  data. 
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of  homogeneity.   As  a  consequence,  the  degree  of  freedom  for 
the  composite  spectrum  is  greatly  increased,  which  implies  a 
confidence  higher  than  95%*      It  is  still  possible  that  this 
periodicity  may  be  due  to  some  procedure  inherent  to  the  NMC 
objective  analysis,  such  as  the  use  of  the  12  hour  forecast 
as  the  first  guess  or  the  scheme  of  adjusting  anomalous  reports 
to  climatology.   Due  to  the  lack  of  an  adequate  number  of  con- 
ventional stations  in  the  equatorial  South  Atlantic  region, 
the  possibility  that  this  4  to  6  day  quasi-periodicity  is 
largely  a  result  of  the  NMC  analysis  cannot  be  completely 
ruled  out 0   However,  in  view  of  the  facts  that  perturbations 
with  the  same  periodicity  have  been  established  beyond  doubt 
in  the  Pacific,  Caribbean  and  North-African  regions  and  that 
the  same  periodicity  is  found  at  the  lower  levels  of  the  two 
stations  used  in  this  study,  it  is  felt  that  there  is  a  reason- 
able degree  of  certainty  about  their  existence  in  the  equatorial 
South  Atlantic,   In  any  case,  further  study  would  be  worth- 
while especially  if  better  data  can  become  availablec 
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FIGURE  2  -  Estimate  of  the  response  function  of  the  high- 
pass  filter. 
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3(a)    -  Autocorrelation  function   (full   line)    and  the 
n  RiJ   (daslied  line).      See  text   for  details. 


FREQUENCY  (  CYCLES  /  DAY  ) 

FIGURE  3(b)  -  Variance  spectrum  (knot2  per  2  if  /50  day"1)  of 
the  1000  mb  v  series  at  Ascension  Island.   Dashed  line  is 
the  95$-  confidence  limit. 
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FIGURE  4(a)  -  Variance  spectra  (knot2  per  2  Tf  /50  day  1)  of  700 
mb  grid  point  winds  at  the  equator.   The  95%   confidence  limit 
(dashed  line)  is  plotted  in  the  0.16-0.24  cpd  band  whenever  the 
variance  exceeds  the  limit. 
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FIGURE  4(b)  -  Same  as  FIGURE  4(a),  except  for  700  mb  and  5°S 
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FIGURE  4(c)    -  Same   as  FIGURE  4(a),    except    for  700  mb   and   10°S 
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FIGURE   5(a)    -  Same   as  FIGURE  4(a),    except   for   500  mb   and  the 
equator. 
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FIGURE  5(b)  -  Same  as  FIGURE  4(a),  except  for  500  mb  and  5°S. 
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FIGURE   5(c)    -  Same   as  FIGURE  4(a),    except    for   500  rab   and   10   S 
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FIGURE   6(a)    -  Same  as  FIGURE  4(a),    except   for   200  mb   and  the 
equator. 
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FIGURE  6(b)    -  Same   as  FIGURE  4(a),    except   for  200  mb   and   5°S. 
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FIGURE   6(c)    -  Same   as  FIGURE  4(a),    except    for   200  mb   and   10°S. 
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FIGURE  7(a)  -  Composite  spectra  (knot2  per  2  TT  /50  day"1)  for 
the  equator.   The  95$  confidence  limit  in  the  band  0.16-0.24 
cpd  is  drawn  v/henever  the  variance  in  the  band  exceeds  the 
limit . 
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FIGURE  7(b)  -  Same  as  FIGURE  7(a),  except  for  5  S. 
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FIGURE  7(c)  -  Same  as  FIGURE  7(a),  except  for  10°S. 
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FIGURE  8  -  Vertically  averaged  spectrum  (m  sec  '  per  2/50  day   ) 
for  the  12  lowest  levels  of  the  v  series  at  Parnaiba  Airport. 
The  95$  confidence  limit  is  depicted  by  the  dashed  line. 
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FIGURE  10  -  Inter-level  phase  relationship  at  Parnaiba 
Airport.   The  base  series  is  801  m.   All  the  corresponding 
coherence  square  values  meet  the  95$  confidence  requirement. 
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FIGURE  11(a)  -  Inter-longitude  phase  relationships  at  the 
equator.  The  base  series  is  15°W.   Positive  phase  differences 
indicate  that  the  base  series  leads  the  others.   Only  values 
with  corresponding  coherence  square  exceeding  the  95> 
confidence  limit  are  plottedo 
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FIGURE  11(b)  -  Same  as  FIGURE  11(a),  except  for  5°S. 
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FIGURE  11(c)  -  Same  as  FIGURE  11(a),  except  for  10°S 
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